function z=rosenbrock(x1,x2)
  z=100*(x2-x1.^2).^2+(1-x1).^2;
end